Mechanisms for the superconducting state from a one-particle derivation of the BCS 

gap equations. 
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The BCS results for the superconducting gap A and Tc are obtained from a one-particle model. 
Superconductivity appears when the electronic energy gains of the band structure surpass the energy 
needed for atomic vibrations or magnetic moment oscillations. The vibration/oscillation amplitudes 
determine the superconducting gap, and the Fermi surface is important for the q-dependence. This 
permits for complementary interpretations of the parameters for superconductivity and modeling of 
density-of-state effects. It also makes the superconducting mechanism less exotic. 

PACS numbers: 74.20.-z,74.20.Fg,74.20.Pq 
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I. INTRODUCTION. 

The Bardeen-Cooper-Schrieffer (BCS) theory is the ba- 
sis for a microscopic understanding of superconductiv- 
ity It has, together with band theory, been success- 
ful for an understanding of the correlation between the 
superconducting (SC) Tc and material parameters like 
phonon frequency, electron-phonon coupling, A, and the 
electronic density-of-states (DOS) at the Fermi energy, 
N(E F ) 0, i]. The good understanding of the " gap ra- 
tio" , the ratio between the SC gap, A, at T—0 and ksTc, 
is another example of the success of the BCS theory. The 
observation of isotope effects are often in good agreement 
with the predictions from BCS, which confirms the de- 
pendence on lattice stiffness. Although the variations of 
Tc with pressure can be complex they are usually under- 
stood from lattice hardening and the behavior of the ma- 
trix elements for electron-phonon coupling [4]. However, 
there are complications Q ; Competition with magnetism 
and electron-electron correlation is reducing Tc- High- 
Tc cuprates are not understood, since their A's are not 
very large. The isotope effect is sometimes very complex, 
even in elementary metals. Other shortcomings might be 
that the many-body BCS formula does not allow for an 
intuitive comprehension of the superconducting mecha- 
nism and there is no direct explanation of the Meissner 
effect [|. 

Here is presented a simple one-electron like deriva- 
tion of the BCS gap equation. Many assumptions are 
the same as in BCS, but some interpretations are differ- 
ent. The many-body formulation of BCS gap equation is 
based on exchange of virtual phonons and the SC state 
appears "when the energy difference between the elec- 
trons involved is less than the phonon energy" In 
the present one-particle derivation the gain in electronic 
energy is caused by dynamic changes of the band struc- 
ture. A one-particle band will have a gap because of a 
periodic potential perturbation, as in the appearance of 
a gap for semiconductors A spontaneous lattice vi- 
bration creates the perturbation and the electronic states 
at k and -k are affected equally (but there is no exchange 
of phonons). Also states on nearby k-points are involved 
as can be found from the band dispersion. Only phonons 



which modifies electronic states around the Fermi surface 
(FS) are of interest, since changes of the bands far from 
the Fermi energy [Ep) make no change in the total ki- 
netic energy. This puts a constraint on the q- vectors of 
the phonons. The difference between the present deriva- 
tion and original BCS is that the energy gain is calculated 
from all one-particle states near Ep, and the SC state ap- 
pears as soon as this gain is larger than the vibrational 
energy. This allows for further insight to the mechanism 
of SC pairing, since the parameters can be understood in 
a simple way from the band structure. 



II. ONE-ELECTRON MODEL AT T=0. 

Phonons and spin fluctuations in the normal (non-SC) 
state are excited thermally following the thermal Bose- 
Einstein occupation, g(T, cu), of the phonon- or spin wave 
DOS, F(u>) or F m (oj), respectively. The averaged atomic 
displacement amplitude for phonons, u, can be calculated 
as function of T [y,0|- Approximate results for Debye like 
spectra make u 2 — » "Sftwr)/2K at low T and iksT/K at 
high T, where K is an effective force constant, which can 
be calculated from K — Maw 2 — d 2 E/du 2 (E is the 
total energy), where Ma is an atomic mass. The cor- 
responding relations for averaged fluctuation amplitudes 
of the magnetic moments, m, are the same, but without 
the polarization factor 3 and with replacement of K with 
K m = d 2 E/dm 2 , which is constant for harmonic depen- 
dence of the total energy as function of m Thus, for 
the normal state, one can estimate u and m at any given 
T from these relations, if the 'force' factors K and K m 
can be calculated or are known from experiment. 

Phonons and spin waves have an influence on the elec- 
tronic state and its DOS [9(, and may cause a pseudo gap 
close to Ep in the normal state of high-Tp copper oxides 
[Iol |. Phonons make a periodic potential perturbation 
along a chain of atoms, 
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if the phonon propagation is along x with wave vector 
A spin wave makes an analogous perturbation 
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within the spin polarized part of the potential. The po- 
tentials for opposite spins are the same except for a phase 
shift of 7r. The result is an anti- ferromagnetic (AFM) spin 
configuration with wave length given by 2n/q. Phonons 
and spin waves can be considered separately, but sev- 
eral works have shown that an important spin-phonon 
coupling (SPC) exists in the cuprates, which also can 
explain many of their unusual properties [ill S3, EH 13 • 

The following development is based on phonon exci- 
tations, but later it will be seen how things will change 
with spin waves and SPC. 

At very low T there are practically no occupation 
of phonons. The Fermi-Dirac occupation /(e,T) = 
l/(e £ / fcfiT -I- 1) is essentially a step function at Ep (here 
Ep is put at zero). The simple nearly free-electron 
(NFE) model, with the periodic potential, leads to a 
gap Eg = 2Vq in the free electron band at a new "zone- 
boundary" k = G/2 and the general band disper- 
sion as function of k is 

e = \{k 2 + (k Gf ± ^{k 2 - (k - G) 2 ) 2 + 4K) 2 ) (2) 

If the gap appears at Ep for this particular value of k 
there is a gain in kinetic energy, E. 

For band energies near the gap it is convenient to ex- 
press the k-dispersion in terms of a linear e measured 
from Ep — 0, so that e = const. ■ k, where k is measured 
from the zone boundary at G/2. The approximation of 
a linear e as function of k is valid for e << W, the band 
width from the bottom of the band to Ep. Thus, 

s = ±^(e 2 + V 2 (3) 

The normal free electron dispersion, e = e, is recovered 
for Vq = 0, and N(e) = N/\de/de\, becomes constant and 
equal to N. The gapped N is zero for energies within ±Vq 
around Ep. 

The electron-phonon coupling A will enhance the elec- 
tronic DOS for energies ±htu around Ep, where it can 
be written NM 2 /K Q. The matrix element M is zero 
outside the interval ±hu>. Inside the interval it can be 
evaluated as (^*(E F , r)dV(r)/du^(E F , r)), which is the 
first order change in energy caused by the perturbation 
dV(r) for du — > 0. For a finite value of u the change in 
energy will be finite and equal to the gap Vq, since Vq/u 
is constant for harmonic vibrations. Thus, instead of cal- 
culating M as a matrix element it is possible to take the 
value directly from the band gap, and M can be written 
Vq/u for energies close to Ep. The separate values of Vq 
and u are important variables, and later it can be verified 
that the band gap Vq is linked to the SC gap A. 

The energy of an atomic oscillation consists of elas- 
tic and kinetic contributions, but its time dependence, 
2U(t) — Ku 2 cos 2 {ijjt) + Maui 2 sin 2 (ujt), is a constant in 
the harmonic approximation. For t — all the energy 
is of elastic origin so the cost in energy to generate an 
atomic vibration can be written U = ^Ku 2 , where u 



refers to the maximal atomic displacement. Totally there 
will be a gain in energy if \U\ < \E\. The system will 
spontaneously generate phonons in such a case, and this 
is the condition for the SC state. Both U and E are the 
energies per unit cell for which the DOS of the normal 
state is N. The condition \U\ = \E\ is written 

l -Ku 2 = [ e(N(e) ~ N(e))de (4) 

for T = 0, where N(e) is the DOS with the gap and 
N(e), the DOS of the normal state, is assumed constant 
within fku around Ep. The integration is only to fuv 
since A, which will appear later in the equation, is zero 
for energies larger than ±huj. 

With a subst itution e 2 = e 2 + A 2 we obtain TV = 
Nlel/^e 2 - Vq 2 and, 



-Ku 2 = Nede- Ne 2 /J{e 2 - V 2 )de (5) 
2 J-huj J -hu v 

With A replacing Vq this gives 

Ku 2 = N{fuv) 2 + NA 2 ln(2fuv/A) - N(hw) 2 (6) 

and 

A = 2hue-^ x (7) 

since A = NA 2 /Ku 2 when the gap at T = is A. 

Therefore, on one hand it can be argued that Vq has 
to be equal to A, the SC gap, since this derivation then 
reproduces the BCS result. But the equivalence between 
A and Vq can also be understood from the fact that Vq is 
a measurable band gap in the superconductor. Further, 
a constant A implies that the phonon amplitude u is pro- 
portional to A, so that u is largest at T=0, and u — > 
when A -> at T -> T c . 

The integral and the interpretations look a bit different 
from BCS but the final result for A is the same. 

Any system will generate phonons as soon as the gain 
in electronic energy generated by the phonon is larger 
than the energy needed for the phonon. In reality other 
more subtle effects, electron-electron correlation energy 
and potential terms will be added to the energy costs and 
those terms will prevent a SC gap in many systems. 

A numerical solution needs some care for the diverging 
part of N, but models with non-constant N and with 
larger ratios of can be studied. 

III. THE LIMIT A ^0. 

The model for finding Tc is obtained from eq. 01 but 
with the Fermi-Dirac function as the T-dependent weight 
factor for N(e) and N(e), and with the integration in 



3 



the interval [—fvjj,tkJ\. Thus, at Tc it is required that 
the phonon energy is equal to the difference between the 
kinetic energy of the gapped and the normal electronic 
DOS, but with the constraint that A — > 0. This is solved 
numerically from 

Ku 2 » NI(hoj, A, T) (8) 

or 

1/A«/(^,A,T)/A 2 (9) 

where 

/hco fhu)' 
efde- / e |e|/V(e 2 - A 2 )fde (10) 
-hw J —huj' 

for A — > 0. The ' in the first integral means that the 
energies where |e| < A are excluded. A factor of A 2 is 
extracted from the right hand side of eq. [8] and combined 
with N and Ku 2 on the left hand side to give A as before. 

The original BCS expression for Tc, which is derived 
from 

r huJ de 

X/gN= / —tanh(e/2k B T) (11) 
Jo e 

(where gN is the coupling constant [151 ]) is independent 
of A. Here, we solve this equation numerically with the 
same precision as the one-particle expression for very 
small A. The results, shown in Fig. 1, tend towards the 
analytic solution of eq. QTJ the well-known BCS formula, 
which can be written 

k B T c = l.l3hue- 1/x (12) 

when gN — A. For example, it can be verified from 
this formula and Fig Q] that a A of 0.5 makes k B Tc ~ 15 
meV when ftio is 100 meV. 

A non-constant DOS is considered by solving eq. [9lfTUl 
with the N(e) directly from eq. 2 with tiw being a large 
fraction of W. The requirement of constant number of 
electrons as function of T in the normal and SC state 
makes the numerical stability more difficult. But the 
results indicate that Tc is lower if the non-constant (free 
electron like) DOS is included in the integrations over 
±Huj compared to the result with a constant DOS. 

For AFM spin waves there is a cost in magnetic en- 
ergy, which in the harmonic approximation can be writ- 
ten U m = ^K m m 2 . The change in potential on some 
site, V m , is positive for one spin and negative for the 
other spin direction. This defines a A s / = NV^/ Km 2 
as a coupling constant for spin- fluctuations [8j. The rest 
of the equations are valid with A s / replacing A and with 
hw s f being the energy of the spin wave. 
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FIG. 1: The full line shows 1 /A as function of fesT as obtained 
for eq. 111! i.e. the numerical BCS result for hu=0.1 eV. The 
thin and broken lines show the corresponding result from eqs. 
[9lfT0l for 5 different A's (0.5, 1, 2, 4 and 8 meV). The result 
for the smallest A is indistinguishable from the BCS result 
over this temperature range. 



IV. SPIN-PHONON COUPLING. 

As was mentioned above, typical atomic displacements 
and magnetic moments from phonons and spin waves in 
the normal state can be determined from thermal excita- 
tions via the effective force constants K and K m . Mag- 
netic moments, with a tendency for a pseudogap in the 
DOS, are driven by thermal excitations, but the left hand 
side of eq. [5] is larger than the right hand side. However, 
the situation might be reversed at lower T, when m is 
supported by the SC gap. In the SC state u and m 
are proportional to the SC gap; u = y/(N/K\)A and 
711 = \/~(N/K m \ s f)A, respectively. Increasing ampli- 
tudes of u in superconductors at low T should be mea- 
surable, but their values are small in conventional super- 
conductors. For example, in Nb, with A f» 1.2, N w 
OJieVatomspin)- 1 , K w 6eV/A 2 Q and A w 3meV, 
u/aQ (ao the lattice constant) will be less than one order 
of magnitude smaller than u from zero-point motion. The 
complex FS of d-band superconductors implies a multi- 
tude of active q- vectors. Lithium, which under pressure 
can have a high Tc and a simple free electron FS, should 
be more promising for detection of u for a few q- vectors. 

Copper oxides with high Tc and relatively simple 
2-dimensional FS, should have sizable amplitudes of 
the magnetic moments (assuming that spin fluctuations 
are responsible for superconductivity), and simple q- 
dependence. Both the superconducting gap and the pseu- 
dogap for T > Tc are sensitive to spin waves (or phonons) 
with the same g-vectors. For instance, fluctuations in 
form of spin waves and SPC in the cuprates are ther- 
mally excited at large T and contribute to a pseudogap 
for T <T* [12]. The fact that m (or u) is proportional to 
A shows that these magnetic fluctuations will reappear 
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in the SC state and become stronger as T — > 0. This is 
in line with the observations of increasing peak intensity 
of spin waves at or below Tq in experiments of inelastic 
neutron scattering on underdoped YBCO 17]. Theoret- 
ical estimates of m, in the range 0.1-0.2 /j,b/Cu in the 
SC state and in the normal state at large T, are not very 
precise because of uncertainties in density functional cal- 
culations. However, the q-dependence is not expected to 
change as T goes below Tc, so the results for spin exci- 
tations calculated for the normal state in ref. [l8| can be 
carried over to the SC state. 

The standard propositions for higher Tc is to increase 
u> (through isotope shifts) and/or A (through higher abso- 
lute value of the DOS [IJI). It is seen that A can remain 
constant and lead to larger A if m (or u) is increased. 
This might be achieved through anharmonicity so the 
maximum m are increased without large changes of K m . 
Another possibility is to modify the way ksT is occupy- 
ing and depleting states around the SC gap and around 
Ep in the normal state. Temperature and the Fermi- 
Dirac distribution cannot change, but different energy 
variations of the DOS near Ep are possible. The effect 
of a gap originating from a potential perturbation on the 
real band structure has to be determined and inserted in 
eqs. IWTUlfor sorting out the effects of non-constant DOS, 
N(e). 

The case with strong SPC in the cuprates leads to large 
enhancements of A s / for some phonons with particular 
g-vectors. Anharmonicity is also expected in SPC when 
the magnetic moments can be enhanced at large u, which 
then lead to mutual softening of phonons and spin waves. 
This can be taken into account in quasi-harmonic vibra- 
tions as amplitude dependent K m and K. When the 
gap at Ep is caused by a few periodicities of phonons or 
spin waves it might be important to consider decoupled 
g-dependencies in the numerical search of A and Tq- For 
instance, large Xspc for selective q on the 2-dimcnsional 
FS of the Iiigh-T£7 cuprates is favorable to d-wave pairing, 
because the strength of the gap is different at different 
parts of the FS |ll|, |l6|. However, the present derivation 
does not yet consider the phase of the gap. In direct 
SPC, the excitations of a phonon and a spin wave are 
made together, which can make u) s f ~ lu and a large 
Aspc- Indirect, weaker X s f are possible for independent 
spin fluctuations which are enhanced by the presence of 
normal thermally excited phonons at large T [l^, [2p| ] . 
This relative strength of SPC is consistent with effective 
couplings derived from optical spectra showing that the 
strongest couplings appear near phonon energies, while 
weaker couplings exist at higher energies [2lj . 



V. FIELD DEPENDENCE. 

From the discussion about selective q-dependence 
there is a possible reason to why a weak magnetic field is 
expelled in a superconductor. As mentioned, one partic- 
ular phonon, q, is responsible for the gap on the para- 




FIG. 2: The full, broken and semi-broken lines show the val- 
ues of the energy integral D(h, T) (see eq. I15[) as function of 
h for low, intermediate, and high temperature, respectively. 
The values are normalized to 1 for h = 0. These calculations 
are made with A = 5mRy and fkv — 50mRy. 

magnetic FS at Ep. A magnetic field, H, will split 
the FS into two, one for each spin ("up" or "down"), 
Ep d = Ep ± (IbH. If SC gaps should appear optimally 
on these two FS, it requires two independent phonons. 
One produces a potential perturbation at e~ l ^~ s * x , and 
the other at e ~ l i?+ 8 )- x ; where 5 is determined by the band 
dispersion and fisH. The sum of these two potentials is 
2cos(8 ■ x)e~ lq ' x , and therefore, even if there is a modula- 
tion given by the cosine function, the effective q remains 
the same and will not fit to the optimal values for the 
two FS. The resulting SC gaps do not appear at Ep on 
the two spin-split bands. This will reduce the gain in 
energy of the SC state, as can be seen from a two-level 
model: Suppose that the DOS of the SC state consists of 
two peaks at ±A (each containing one electron of each 
spin) and that Ep—0 in zero field. The kinetic energy in 
this case 

£o = 2(/(-A,T)(-A) + /(A,T)(A)) (13) 

A field puts the states asymmetrically around Ep and the 
kinetic energy for the two spin states ("up" or "down") 
will be 

E% = f(-A±h,T)(-A±h) + f(A±h,T)(A±h) (14) 

where h — ^jlbH. The result is that E < E^ + E^ for 
most T > (but for ksT and h being small in comparison 
to A), i.e. the symmetric state with no field has the 
lowest energy. 

The increasing kinetic energy can also be demonstrated 
for the approximation of a constant DOS by adding and 
subtracting the field h in the Fermi-Dirac function. Fig- 
ure shows the energy difference, D(h,T), of kinetic 
energy for the gapped superconducting DOS with and 
without field, which is calculated as: 
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D(h,T)= / eN(e)(f(e+h,T)+f(e-h,T)-2f(e,T))de 
J — huj' 

(15) 

when h < A for low and high T (« A). The increase of 
D as function of the field h is because the thermal occu- 
pation can be made more efficiently if Ep is closer to the 
DOS peak (on N) above the gap for "majority" and closer 
to the DOS peak below the gap in the "minority" states, 
than if Ep is in the middle of the gap. Thus, the gapped 
state with /xsiJ=0 has the lowest kinetic energy. This 
state will be preferred by the system as long as screening 
of an external held can be made through superconducting 
currents. The model shows that the minimum at h = 
is less profound for large T, or when fcgT exceeds about 
5A. The feedback from the transfer of minority to ma- 
jority spin states and effects of a non-constant DOS are 
not included in the model. 



VI. CONCLUSION. 



The BCS formulas for A at T = and T c at A = 
are derived directly from the one-particle DOS functions 
of the gapped and normal state band structures. This 
allows for an easy comprehension and further interpreta- 
tions of the SC mechanism. While phonons and/or spin 
waves are excited thermally in the normal state, they are 
generated via the electronic band gap in the supercon- 
ducting state. Atomic displacements of harmonic vibra- 
tions and magnetic moments of harmonic spin fluctua- 
tions are proportional to the SC gap. Since the SC gap 
is closely related to the gap of the perturbed band struc- 
ture, it will be interesting to consider DOS functions in 
materials with impurities and other defects via supercell 
calculations. 
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